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Abstract 

The spherically symmetric Einstein- Vlasov system is considered in 
Schwarzschild coordinates and in maximal-isotropic coordinates. An 
open problem is the issue of global existence for initial data without 
size restrictions. The main purpose of the present work is to propose 
a method of approach for general initial data, which improves the reg- 
ularity of the terms that need to be estimated compared to previous 
methods. We prove that global existence holds outside the centre in 
both these coordinate systems. In the Schwarzschild case we improve 
the bound on the momentum support obtained in |19j for compact ini- 
tial data. The improvement implies that we can admit non-compact 
data with both ingoing and outgoing matter. This extends one of the 
results in j6]. In particular our method avoids the difficult task of 
treating the pointwise matter terms. Furthermore, we show that sin- 
gularities never form in Schwarzschild time for ingoing matter as long 
as 3m < r. This removes an additional assumption made in [1 . Our 
result in maximal-isotropic coordinates is analogous to the result in 
[20] , but our method is different and it improves the regularity of the 
terms that need to be estimated for proving global existence in general. 

1 Introduction 

In the present work we investigate the issue of global existence for the spheri- 
cally symmetric Einstein- Vlasov system when the initial data is unrestricted 
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in size. The system is analyzed in Schwarzschild coordinates, i.e. in a po- 
lar time gauge, and in maximal-isotropic coordinates where a maximal time 
gauge is imposed. These coordinate systems are often, in the literature, con- 
jectured to be singularity avoiding. However, there is to our knowledge no 
proof of this statement for any matter model and it would be very satisfy- 
ing to provide an answer to this conjecture for the Einstein- Vlasov system. 
Moreover, a proof of global existence would be of great importance due to 
its relation to the weak cosmic censorship conjecture, cf. [9[ IT U [13 ], [Tj 14"] 15]. 
A third motivation for our interest in this problem is the fact that global 
existence for general data have been obtained for the Vlasov-Poisson system, 
i.e. Newtonian gravity coupled to Vlasov matter. Batt [8] showed global 
existence in the spherically symmetric case 1977, and the general case was 
settled independently by Pfaffelmoser [16], and Lions and Perthame [14] in 
1991. It is thus natural to ask if similar results can be obtained when New- 
tonian gravity is replaced by general relativity. One should of course bear 
in mind that the situations are fundamentally different since in the latter 
case there exist data which lead to spacetime singularities, cf. [H [2Tj . 
Nevertheless, as mentioned above, global existence may hold for polar or 
maximal time slicing. 

The issue of global existence for the spherically symmetric Einstein- 
Vlasov system has previously been investigated in several papers, cf. 
[HIHESCE1II21E1I11E]. Global existence for small initial data has been 
proved in [18] and [12] for massive and massless particles respectively. In [1] 
initial data which guarantee formation of black holes are constructed, and 
it is proved that for a particular class of such initial data, with a steady 
state in the interior of the surrounding matter, global existence holds in 
Schwarzschild time. In [3j global existence is shown to hold in a maximal 
time gauge for rapidly outgoing matter. The methods of proofs in these 
cases are all tailored to treat special initial data and they will likely not 
apply in more general situations. The investigations [191 EH CO E] are con- 
ditional in the sense that assumptions are made on the solutions, and not 
only on the initial data. However, the methods are general and cover large 
classes of initial data and can thus be thought of as possible approaches 
for treating the general case. The main purpose of the present work is to 
propose a method of approach for general initial data, which improves the 
regularity of the terms that need to be estimated compared to the methods 
in |19 [ 1201 \T\ l6j. Furthermore, the method improves and simplifies some of 
the previous results. To be more precise, let us discuss the relation between 
these studies. 

An important quantity in the study of the Einstein- Vlasov system and 
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the Vlasov-Poisson system is 

Q{t) : = sup{l + \v\ : 3(s,x) € [0,t] x M 3 such that f{s,x,v) ^ 0}. (1.1) 

Here / is the density function on phase, and u 6 M 3 is the momentum. Q(t) 
measures the support of the momenta, and the content of the continuation 
criterion for these systems, cf. [13 [201 [8], 1S that solutions can be extended 
as long as Q(t) remains bounded. The definition of Q only applies in the case 
when the initial data have compact support in the momentum variables. For 
non-compact initial data, it was shown in [6j that solutions can be extended 
as long as Q(t) is bounded where 

Q(t) :=sup| 1 + l ^ ( ^°| r ^ )l \ 0<s<t, M)€su P p/l, (1.2) 

Here V is a solution of the characteristic system, cf. (|2.18p . The analysis 
in |19} IU [6] is carried out in Schwarzschild coordinates. The main result in 
|19j shows that as long as there is no matter in the ball 

{iGR 3 : \x\ < e}, 

the estimate 

Q(t) < e l°gQ(0)e^ (L3) 

holds. Here C(e) is a constant which depends on e. In view of the continu- 
ation criterion this can thus be viewed as a global existence result outside 
the centre of symmetry for initial data with compact support. The bound is 
obtained by estimating each term individually in the characteristic equation 
(I2.18P for the radial momentum. This involves a particular difficulty. Let us 
consider the term involving p r in (12.18p . The Einstein equations imply that 

Mr = ^e 2A + 4vrrpe 2A =: T x + T 2 , 

where m is the quasi local mass and p is the pressure. There is a distinct 
difference between the terms T\ and T 2 due to the fact that m can be 
regarded as an average, since it is given as a space integral of the energy 
density p, whereas p is a pointwise term. Also the term involving in (|2.18p 
is a pointwise term in this sense. The method in [19] is able to estimate the 
pointwise terms outside the centre but generally it seems very unpleasant 
to have to treat these terms. In the present work we give an alternative 
and simplified proof of the result in [19J. In particular our method avoids 
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the pointwise terms by using the fact that the characteristic system can be 
written in a form such that Green's formula in the plane can be applied. This 
results in a combination of terms involving second order derivatives which 
can be substituted for by one of the Einstein equations. This method was 
first introduced in pQ but here the set up is different and the application 
of Green's formula becomes very natural. In addition the bound of Q is 
improved compared to (jl.3jl and reads 

Q(t)<(Q(0) + ^)e c ^ 1+t ^. 

This bound implies that also Q(t) is bounded and therefore a consequence 
of the method is that global existence outside the centre also holds for non- 
compact initial data. This improves the result in [6] where only non-compact 
data for which the matter is ingoing, and such that the matter keeps on going 
inwards for all times, are admitted. Another consequence of our method is 
that global existence holds for ingoing matter as long as 3m(t, r) < r. Note 
that in Schwarzschild coordinates 2m(t, r) < r always, and that there are 
closed null geodesies if 3m = r in the Schwarzschild spacetime. This result 
was already proved in pQ but an additional assumption was imposed which 
now has been dropped. 

We now turn to the case of maximal- isotropic coordinates. Rendall shows 
in [20] global existence outside the centre in maximal-isotropic coordinates. 
The bound on Q(t) is again obtained by estimating each term in the char- 
acteristic equation, cf. equation (16.171) . In this case there are no pointwise 
terms in contrast to the case with Schwarzschild coordinates. The terms are 
however, in analogy with the Schwarzschild case, strongly singular at the 
centre. The method that we use in Schwarzschild coordinates also applies 
in this case. The improvement lies in the fact that there is some gain in reg- 
ularity, i.e. the terms that need to be estimated are less singular compared 
to the corresponding terms in |20j . Indeed, these terms can schematically, 
in both the Schwarzschild and maximal-isotropic case, be written as 

/ g(t,v)dv, (1.4) 
Jo 

and the known a priori bounds read 

vg(t,v) d v- (1-5) 

Hence, the degree of the singularity is of order one in the radial variable. 
Roughly, for the methods in [19] and [20] the singularity is of second order. 
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An advantage of our method is clearly that it applies in both cases, and it 
also turns out that the principal term to be estimated in the two cases is the 
spacetime integral of the Gauss curvature of the two dimensional quotient 
manifold M/SO{3), where Ai is the four dimensional spacetime manifold. 
This is interesting since the Gauss curvature is a coordinate independent 
quantity. Moreover, the principal term in |10j . p. 1172, has the same form. 
Finally we mention that our method also applies in the case of maximal- 
areal coordinates which were used in the study [3] . The analysis in this case 
is completely analogous and is left out. 

The outline of the paper is as follows. The Einstein- Vlasov system in 
Schwarzschild coordinates is treated in sections 2-5. In section 2 the system 
is formulated and the a priori bounds are given in section 3. Section 4 is 
devoted to the proof of global existence outside the centre and in section 
5 global existence is shown for ingoing matter which satisfies m/r < 1/3. 
Sections 6-8 concern the system in maximal-isotropic coordinates. The sys- 
tem is given in section 6, and in section 7 the necessary a priori bounds are 
derived. Section 8 is devoted to the proof of the global existence theorem 
outside the centre. 



2 The Einstein-Vlasov system 

For an introduction to the Einstein-Vlasov system and kinetic theory we 
refer to [21 [20], and for a careful derivation of the system given below we 
refer to |17] , In Schwarzschild coordinates the spherically symmetric metric 
takes the form 

ds 2 = _ e W,r) dt l + e m,r) dr 2 + + ^2 ^ 

The Einstein equations read 

e~ 2A (2rA r - 1) + 1 = 8irr 2 p, (2.2) 
e- 2X (2rfi r + 1) - 1 = 8vrr 2 p, (2.3) 
\ t = -4irre x+ vj, (2.4) 

e- 2X {n rr + (ji, - X r ){p r + -)) - e~ 2 ^{\ tt + A t (A t - ih)) = 8tt Pt . (2.5) 

r 

The indicies t and r denote derivatives. The Vlasov equation for the density 
distribution function / = f(t,r,w,L) is given by 

d t f + e»- x ^-d r f - (X t w + e^ X fM r E - e ^ x ^-)d w f = 0, (2.6) 



5 



where 

E = E(r, w, L) = + + L/r 2 . (2.7) 

Here w 6 (—00, 00) can be thought of as the radial component of the mo- 
mentum variables, and L £ [0, 00) is the square of the angular momentum. 
The matter quantities are defined by 

POO /"OO 



r-co r-co 

p(t, r ) = —l I Ef(t,r,w,L) dwdL, (2.8) 
r J-00 Jo 

roc rco ^2 

P^ r ) = ^i / - W f(t,r,w,L)dwdL, (2.9) 

J — oo J 
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-oo JO 

/OO /'OO 
/ wf(t,r,w,L), dwdL, (2.10) 
-oo JO 
^ roc roo j 

PT(t,r) = — [ / -f(t,r,w,L)dwdL. (2.11) 

zr J-oo JO & 

Here p,p,j and pt are the energy density, the radial pressure, the current 
and the tangential pressure respectively. The following boundary conditions 
are imposed to ensure asymptotic flatness 

lim X(t,r) = lim fj,(t,r) = 0, (2.12) 

r— >oo r— Yoa 

and a regular centre requires 

A(t,0) = 0, t > 0. (2.13) 

We point out that the Einstein equations are not independent and that e.g. 
the equations fl21J and ([23]) follow by $T2§, §23$) and (|23|) . 

o 

As initial data it is sufficient to prescribe a distribution function / = 
f(r, w,L) > such that 

Wp( V )d V <^. (2.14) 
Here we denote by p the energy density induced by the initial distribution 

o 

function /. This condition ensures that no trapped surfaces are present 

O 

initially. Given /, equations (|2.2j) and (|2.3p can be solved to give A and p 

o 

at t = 0. We will only consider initial data such that / = if L > L2, for 

o 

some L2 > 0, and such that /= if r > R2, for some R2 > 0. If in addition 
the initial data is C 1 we say that it is regular. In most of the previous 
investigations the condition of compact support on the momentum variable 
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w has been included in the definition of regular data. The exception is [6] 
where non-compact initial data is studied and a decay condition replaces 
the assumption of compact support. In this study we also include data with 
non-compact support and we impose the decay condition from |6] 

sup \w\ f(r, w,L) <oo. (2-15) 

(r,to,L)6R 3 

We distinguish between the two cases; regular initial data with compact 
support and regular initial data which satisfy the decay condition (j2.15|) . 
We denote these classes of initial data by I c and TP respectively. 

The main results below concern subclasses of T c and TP . Given R\ > 0, 
we define the subclass T C (R\) with radial cut-off by 

X C (R l ) = {/£ X C : /= for r < R x }. 

The subclass l D (Ri) is defined analogously. 

Let us now write down a couple of facts about the system (|2.2p - (|2.13p . 
A solution to the Vlasov equation can be written 

f(t, r, w, L) = f (R(Q, t, r, w, L), W(0, t, r, w, L),L), (2.16) 

where R and W are solutions of the characteristic system 

dR _ r (M-A)(s,fl) W . , 

ds E(R,W,LY [ ' 

^- = -\ t (s,R)W -e^- x ^ s ' R) fi r (s,R)E(R,W,L) 
ds 

+e R3E(R,W,L)> [2A *> 

such that (R(s, t, r, w, L), W(s, t, r, w, L), L) = (r,w,L) when s = t. This 
representation shows that / is nonnegative for all t > 0, ||/||oo = ||/o||oo, 
and that f(t, r, w, L) = if L > L2. The quasi local mass m is defined by 



m(t, r) = 4tt / ^p(t,rj)drj, (2.19) 
J 

and by integrating (|2.2p we find 
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A fact that we will need is that 



H + A < 0. 

This is easily seen by adding the equations (|2.2p and (j2.3|) . which gives 



and then using the boundary conditions on A and [i. Furthermore, from 
(|2,20p we get that A > 0, and it follows that \i < 0. We also introduce the 
notations ft and jl. From equation f|2.3|) we have 



Finally, we note that in [18] and [B] local existence theorems are proved for 
compact and non-compact initial data respectively, and it will be used below 
that solutions exist on some time interval [0, T[. 

3 A priori bounds in Schwarzschild coordinates 

In this section we collect the a priori bounds that we need in the proofs 
below. There are two known conserved quantities for the Einstein- Vlasov 
system, the number of particles and the ADM mass M. Here, we will only 
need the latter which is given by 



The conservation of the ADM mass follows from general arguments but it 
can easily be obtained by simply taking the time derivative of the integral 
expression and use of the Vlasov equation and the Einstein equations. The 
following results are given in pQ but since the proofs are very short we have 
included them here for completeness. By a regular solution we mean a 
solution which is launched by regular initial data. 

Lemma 1 Let (/, [A, A) be a regular solution to the Einstein- Vlasov system. 



X r + [A r > 




(2.21) 




(3.1) 



Then 




(3.2) 



(3.3) 



S 



Proof. Using the boundary condition (|2.12p we get 

r 00 d 

1 >l- e » +x (t,0) = / — e^ +x dr 
Jo dr 

(jir + \ r )e» +x dr. 

The right hand side equals (|3,2p by equations (|2.2p and (|2.3p which completes 
the first part of the lemma. The second part follows by studying instead 
of e" +A . 

□ 

Next we show that not only p(t, •) £ L , which follows from the conservation 
of the ADM mass, but that also e 2X p(t, •) G L 1 . 

Lemma 2 Let (/, /x, A) 6e a regular solution to the Einstein- Vlasov system. 
Then 

poo roc ^ 

/ r 2 e 2X p(t,r)dr < / r 2 e 2X p(0,r)dr -\ . (3.4) 

Jo Jo 87r 

Proof. Using the Vlasov equation we obtain 

d t (r 2 e 2X p(t,r)) = -d r (r 2 e» +x j) - re» +x 2je 2X ™ 

< -d r (r 2 e^ +A j) + Ke tx+x {p + p)e 2X . 

Here we used that m/r < 1/2 together with the elementary inequality 2|j'| < 
p + p, which follows from the expressions (|2.8j) - (|2.10j) . In view of (|2.16p and 
()2.17p we see that limj._5.00 r 2 j(t, r) = 0, since the initial data has compact 
support in r. Since the solution is regular and hence bounded the boundary 
term at r = also vanishes (as a matter of fact spherical symmetry even 
implies that j(t, 0) = 0). Thus, by lemma 1 we get 

d f°° 1 
which completes the proof of the lemma. 

□ 
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4 A regularity result in Schwarzschild coordinates 



Theorem 1 Let < e < R\. Consider a solution of the spherically sym- 
metric Einstein- Vlasov system, launched by initial data in X on its 
maximal time interval [0, T[ of existence. If /(s,r, •,•) = 0, for (s,r) £ 
[0,i[x[0,e], then 

Ct ca+t) 

Q(t) < (1 + -)e— . (4.1) 
e z 

In particular, if f(t, r, •, •) = for (t, r) G [0, T[x [0, e], i/ien T = oo. 

The last statement in the theorem is a consequence of the continuation 
criterion derived in [5J. The theorem thus improves the result in [B] which 
was restricted to a special class of initial data where all the matter is ingoing. 
The result holds in particular for compactly supported data, and the bound 
(|4.ip improves the bound in [19J, cf. inequality fjl .3|) . 

Proof: We consider the quantities G = E + W and H = E — W, which 
satisfy G > 0, H > 0. Along a characteristic (R(s), W(s), L) we have by 
fl2JH) and (l2J8j) . 



and 



ds 

dH 

ds 



A t — + ^ r e"- A 



G + 



Le^ 



-x 



H 



R 3 E 
Le^~ x 

~pJe~' 



Let us first consider the quantity H. From (|4.3p we have 



dH 

ds 



< 



A t — +yu r e 



//—A 



if. 



It follows that 

ff(i) < J ff(o) e / o[ At(s ' K(s)) ^y + ^ (s ' i?(s))e( ^ A)(s ' fl(s)) ] ds . 

Let us denote the curve (s, R(s)), < s < t, by 7. By using that 



(4.2) 
(4.3) 

(4.4) 

(4.5) 



dR 



W 



dt E 

the integral above can be written as the curve integral 

A-n+W,r) Xt (t,r)dr + e^- x)(t ' r) fi r (t,r)dt. 



(4.6) 
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Let T denote the closed curve r := 7 + Cf + Coo + Co, oriented clockwise, 
where 

C t = {(t,r):R(t)<r<R oc }, (4.7) 
C oo = {(a > /2oo):t>s>0} J (4.8) 
C = {(0,r):R OD >r>R(0)}. (4.9) 

Here R^ > R2 + t, so that / = when r > R^. We now apply Green's 
formula in the plane and use equation (|2.5p to obtain 



y e^ +A A t dr + e^ X fi r dt 



(X tt + (At - fi t ) A t - (//jt + (fi r - A r )^ r )e M A dtdr 
e fi-\ 

(fi r - \ r ) 8irp T e^ +x dtdr 

r 

= j j e» +x - MP-P) ~ torprj dtdr. (4.10) 

Using the hypothesis that matter stays away from the region r < e, we have 
R(s) > e, < s < t, and we get in view of (j3.3[) . together with the facts 
that A > 0, p — p > and q > 0, 



2m \ It 

- 4vr(/9 - p) - &irp T dtdr < — . (4.11) 



r 3 / € 



Remark 1: The integrand above has a geometrical meaning. It is the scalar 
curvature of the quotient manifold Ai/SO(3) with metric ds 2 = —e 2>1 dt 2 + 
e 2X dr 2 . It is interesting to note that the principal term in [10], p. 1172, has 
the same form. 

Now we wish to estimate the curve integral fj4.6f) . We have 

/ ... = (£...- j 4irrje 2X dr- [ j^-dt - [ Aitrje 2X dr. (4.12) 
Jj JT JCt iCoo ^00 JCq 

Here we used that At = 47rrje^ +x and that fi r (t, Roo) = The integral 
along Co is given in terms of the initial data, 

4vrrje 2A dr| < C. (4.13) 

Co 
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By letting i?oo — > oo the integral along Coo vanishes and it remains to 
estimate the contribution from the integral along Ct- In view of (|3.4p we get 

I 1 f'°° C 

1/ <iitrje 2X dr\< - 4Trr 2 p(t, r)e 2X{t ' r) dr < -(1 + t). (4.14) 

Jc t e Jo e 

The esimate of G is very similar. Indeed, by the hypothesis of the theorem 
we have R{t) > e, thus 

-We s - £ V- (4 ' 15) 

We derive in view of (14.21) 



G(t) < G(0)e 
Jo e2 

Note that the integrals in the exponent above are identical to the integral 
in (|4.5p except for the sign. We can accordingly use the same arguments as 
above, the only difference is that we use (|4.14p twice, for Ct and C T , and 
that the integral in ()4.10p has opposite sign. Since 2px < p— p it is sufficient 
to estimate the integral 

8ir(p - p)e fl+x dtdr. (4.16) 

n 

We use the bound (|3.2j) to obtain the estimate ()4.1ip also in this case. Hence, 
we have the following bounds 

Ct Oll+t) 

Gf)<(G0+9eV, (4.17) 

and 

H(t) < H^e 2 ^. (4.18) 
Since E = G + H this implies that 

Ct cn+t) 

Q(t)<(l + —)e—. (4.19) 

e z 

This completes the proof of Theorem [TJ 

□ 
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5 Global existence for ingoing matter with 3m < r 

In this section we consider compactly supported initial data / G X° {R\) such 
that for given < L\ < L2 and P > 0, it holds that L > L±, and w < —P 

o 

for all (w,L) G supp/. Furthermore, the initial data have the property that 
3m < r everywhere. Note that in Schwarzschild spacetime there are closed 
null geodesies when 3m = r. Let us denote this class of initial data by 
l c (Ri, Li,P, 3) C T c {Ri). By continuity there is a T x > such that w < 
for w G supp/(t), and 3m(i,r) < r everywhere, for t < Tj. We will show 
that on the time interval [0,Ti] singularities do not form in the evolution. 
This can be phrased as a global existence result for ingoing matter satisfying 
3m < r. A similar result is proved in [T] but an additional assumption on 
the solution is imposed and the present result is thus an improvement. 

Theorem 2 Consider a solution to the spherically symmetric Einstein- 
Vlasov system launched by initial data f G T°{Ri, L\, P, 3). Let T > 
be the maximal time interval on which the solution exists and let T\ be as 
above. IfT < oo then T\ < T. 

Proof: Let us consider the quantity He" along a given characteristic 
(t, R(t), W(t),L), with W < 0. In view of and (jMTjl we have 



W X 



+ /i* + /^r-^-e' 



He" 



Le"~ 



R 3 E 



(5.1) 



Since fj, = ji + fi we have 



— = /i r .e M ~ A +/i r e /i " A (l+— ) = fi r e"~ x +jl r e^ x 



Since H = E - W and 



me 



2A 



Mr 



£(£ - W) 
(5.2) 



we obtain 
d -(He") = 



ds 

We note that if 



W , 



He -\-e A 5 5— e". 

(5.3) 

(5.4) 



m(t,R(t)) 1 
~ 3' 
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we have 



ds 



m{t,R(t)) 2 \(t,R(t)) = m(t,R(t)) 
R(t) 

m 

{He*) < 



< 1. 



(5.5) 



Hence, as long as (|5.4p holds true we get in view of f|5.3j) the inequality 
d 



W j 
\t— + fi r e^ x + p t 



He* + e^~ x ^-e*. (5.6) 



Since ER > \[L\ the last term is bounded. We can without loss of generality 
assume that He^ > 1 on [0, Ti], and we get for t <T%, 

Xe J?[At(«,fl(.))^+Ar(«,flW)e^- A >(-- H W)]d._ (57) 

Let us denote by 7 the curve (t, R(t)), < t < T\. The last integral in (15.71) 
can be written as 



; (- M +A)(t,r) At ^ r ) dr + e 0*-A)(t,r)^ r ( t) 



(5.8) 



We will apply the Green formula to this curve integral and we introduce as 
above the closed curve T = 7 + Ct x + Coo + Co where Ct x , Cx, and Co are 
defined as in (18.711. We have 



j e-» +x \ t dr + e^ x p r ds 

d t ( e-» +x \ t ) - d r ( e^ x jl r ) dsdr 



d t { e-» +x \ t ) - d r ( e^ x u r ) dsdr + / I d r ( e» +x ^) dsdr. 



(5.9) 



For the first integral above we use (|4.10p and we obtain the identity 

-^ +x X t dr + e^ x jlrds 



— — 4ir(p — p) — 8npT ) dtdr 



+ f f e^ +x f( n r + A r )^ + 4ttp - ^\ dtdr 



4vre^ 



+A 



(p + p)e 2X — + 2p- p 



drds 



+ 




n J -00 Jo 



00 47r 2 eAt +A 
r 2 E 



f(t,r,w,F)dFdwdrds. (5.10) 
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Here we used that \i T + A r = 47rr(p + p)e 2X and that 

poo poo 4 2 /x+A 

2 PT = p-p- / f —f(t,r,w,F)dFdw. (5.11) 

J-oo JO r ^ 

Recall from (|5.7p that the integral involving /i should also be taken into 
account. Since 

POO 

At(t,r)= / 4irje 2X e»+ x dr), (5.12) 

we get 

/ jk(s,R(s))ds = [ [ Anje 2X e^ +x dsdr. (5.13) 
Hence we obtain with Cb '■= Ct x + + Co, 

T 

o 

(p + pje^— + 2p-p + ie 2 
r 



</e ^ +x \ t dr + e^ X fi r ds + [ fi t (s,r)ds 

A 7o 



47^ 



+A 



drds 



P P POO POO ^ 2 n+A 

+ / / / / 5 f(t,r,w,F)dFdwdrds 

J JnJ-ooJo r E 

- [ er fl+x X t dr + e f "- x il r ds. (5.14) 



>C B 

The integral over Cb equals 



p POO 

- e-^ +x X t dr + e^- x fi r ds= 4irrj(Ti, r)e 2X(Tl ' r) dr 

Jc B JRin) 
t t T\ foo 

+ / 47ri2oop(8, R^e^+^'^ds- / 4vrri(0, r)e 2X ^ r Ur. (5.15) 
Jo JR(0) 

The first integral on the right hand side is negative since matter is ingoing, 
the second vanishes when is large since p has compact support and the 
third integral depends only on the initial data. For the second integral in 
(|5.14p we note that 

poo poo a 2 p n+\ 4_ 
/ / -^^f(t,r,w,F)dFdw<—r 2 p, (5.16) 
J-oo Jo r 2 E Li 

which in view of (|3.ip implies that it is bounded. It remains to show that 
the first integral in (|5.14p is bounded. Since 1 + 2^e 2A = e 2X we find that 



15 



the integrand is non-positive 

{p + p)e 2X — + 2p-p + je 2X 
r 

= -{p - p)(l - e 2A ™) + (p + j)e 2A < 0. 

Here we used that p < p, e 2X m/r < 1, and that p + j < in view of f|2.9j) 
and (|2.1(jp since matter is ingoing. Hence f/e^ is bounded on [0,Ti]. Now, 
from the characteristic equation (|2,17p we have that 

±R-i = -^ A < ^=e» <J^L< C( Tl )R-\ (5.17) 
ds R 2 E ~ R^ZJ ~ RJT X ~ 

It follows that for any characteristic (t,R(t),W(t)), R(Ti) > e, for some 
e > 0, and Theorem [1] thus applies which shows that the solution can be 
extended beyond t = T\. Hence, if T < oo, then T\ < T, and the proof of 
Theorem [2] is complete. 

□ 



6 The Einstein-Vlasov system in maximal- 
isotropic coordinates 

In |20] the Einstein-Vlasov system is studied in maximal-isotropic coordi- 
nates where the metric reads 

ds 2 = -a(t, R)dt 2 + A 2 (t, R)[(dR + 0(t, R)dt) 2 + R 2 {d6 2 + sin 2 Odcf)]. 

The condition that the hypersurfaces of constant time are maximal implies, 
cf. [20], that the field equations take the following form: 

dR{R2 WA ] = -I a5/2r2 & 2 + 16 -p) 

2 1 3 

a RR + -^glr + -A R a R = aA 2 (-K 2 + 4vr(p + p + 2p T )) (6.2) 

Kr + 3(jA R + h)K = 8vrj (6.3) 

fa--^P= \ aK ( 6 - 4 ) 

A t = -aKA + d R {pA) (6.5) 
1 1 

K t = --^orr + -^A R a R + /3Kr + A7ta(2p T -p- p) 

+ a[-^A RR + -^{A R f-^A R ] (6.6) 
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The indicies t and R denote derivatives but sometimes we also use dt 
and dp/. In [20J the Vlasov equation is written in different coordinates 
than we use here. The relation is as follows. Let x % be the coordi- 
nates (R sin 9 cos (ft, R sin 9 cos (ft, R cos 9) and define an orthonormal frame 
by &i = A~ 1 d/dx l . In [2DJ the mass shell is coordinatized by [t,x a ,v % ), 
where v l denote the components of a vector in the orthonormal frame, and 
the Vlasov equation in [2D] is accordingly written in these variables. Now, let 
w = (v ■ x) / R, and let L = A 2 R 2 (v 2 — (x ■ v) 2 /R 2 ), then the Vlasov equation 
for the density distribution function / = f(t,r,w,L) takes the form 

dtf + {~~ Wnf + (-^ + *Kw + + ^))d w f = 0, (6.7) 

where 



E- L\R.w.L) - Jl + W 2 + -L-. (6i 



The matter quantities are denned by 

f'OO poo 

P^ R ) = ^pj J o Ef(t,r,w,L) dwdL, (6.9) 



-oo J 
foo roc 



rco roc w 2 

P^ R ) = ^R2j I -^f(t,r,w,L)dwdL, (6.10) 

7T f°° f°° 

j(t,R) = ^2^2 J J wf{t,r,w,L), dwdL, (6.11) 

f'OO pOQ T 

PT{t,R) = 2A*mJ_ J E f{t ^ w ' L)dwdL - (6 ' 12) 

To ensure asymptotical flatness we impose the boundary conditions 

lim A(t,R) = 1, lim a{t,R) = 1, lim p(t,R) = 0, (6.13) 

R— too R—too R—too 

and to ensure a regular centre we require that 

A(t,0) = 1, yt> 0. (6.14) 

In the case of Schwarzschild coordinates the constraint equations are easily 
solved for any given / G I . This is not the case in maximal- isotropic 
coordinates. Here we assume that initial data are given, such that the 
constraint equations are satisfied, and such that < / G C has compact 

o 

support, f(R,w,L) = whenever R < R\, R > R2 or if L > L2, for given 
< R\ < R2 < 00, and L2 > 0. Hence, we assume properties which are 
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similar to the initial data class TP(R\), with the notable difference that the 
condition (|2,14p is not imposed since in these coordinates trapped surfaces 
are admitted. Let us denote this class of initial data by J c {R\). 
For a given metric the solution of the Vlasov equation is given by 

f(t, R, w, L) = f (K(0, t, R, w, L), W(0, t, R, w, L)), (6.15) 

where 1Z and W are solutions of the characteristic system 

dU f aW m rait* 

H = { A^-^ (fU6) 
dW Ea R aL 1 A R 

~~ds = A~ aK W + EAW { K + T>- (6 - 1?) 

Here K(t, t, R, w, L) = R, and W(t, t, R, w, L) = w. 

7 A priori bounds in maximal-isotropic coordi- 
nates 

In this section we collect the bounds needed in the next section. The fol- 
lowing estimates are derived in [20] by using the results of |15j . 

^I<|^<1- (T.1) 



Moreover, it follows from equation (|6.4p . cf. [20J, that 

\/3\ < 3. (7.2) 

Another useful bound follows from conservation of the ADM mass M, which 
is given by, cf. [20] . 

1 f°° 3 
M = - A 5/2 R 2 (-K 2 + 167iy>) dR. (7.3) 
4 Jo 2 

The following estimate plays the role of Lemma 1 in Schwarzschild coordi- 
nates. 

Lemma 3 Let (/, a, A, f3) be a regular solution of the Einstein- Vlasov sys- 
tem /TOD- fO^ ). Then 

/ aA 2 r](-K 2 + 4TT(p + p + 2p T ))dr] < 1. (7.4) 
Jo * 
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Proof: Consider the second order equation (|6.2p for a. The left hand side 
of (|6.2p equals 

' d R {R 2 Aa R ), 



R 2 A 

which implies that 

a R (R) = aA 3 V 2 (^K 2 + 47r(p + p + 2p T ))d7 ] . 

By using the boundary condition a(t, oo) = 1 we get 

l>l_ a (0) = J _( J aiV(/ 2 + 44 P+I) + 2|) T ))^)^. (7.5) 

From equation (|6,ip we get that A is monotonic and decreasing. Changing 
the order of integration and using the monotonicity of A we get 

/"OO /'OO j Q 

l>y (y jp J ^dR)aA 3 r ] 2 (-K 2 + 4Tr(p + p + 2p T ))dr ] 
f°° 3 

= / aA 2 r](-K 2 + 4ir(p + p + 2p T )) d V . (7.6) 

□ 

We also need to establish a decay result for large R. 

Lemma 4 A solution of the Einstein- Vlasov system \6.1\) -( 'Un^ satisfies 

I' 1 a R (s,R) 

/ — — — ds — >■ as R — > oo. (7.7) 

7 A(s,i?) 

Proof: Using (|7.3p together with equation (|6.ip we obtain 

Ar M 

Va~ r 2 ' 

which in view of the boundary condition A(oo) = 1 implies that 

A(t,R)<(l + ^) 2 . (7.8) 

From (|6.16p and (|7.2p we obtain that all characteristics originating from the 
support of the matter have R(t) < R2 + 4i =: R 3 (t). From equation (|6.3I) 
we have in view of (|7.ip 

rR rRs(t) 
\{R 3 A 3 K)(t,R)\ = I / %Txr l 3 A A jdr ] \ = | / 8ttt] 3 A 4 j drj\ 
Jo Jo 

= \(R 3 A 3 K)(t,R 3 (t))\ < 2R 3 (t) 2 A 2 (t,R 3 (t)). (7.9) 
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From (17.81) we have 



Thus, since A > 1, 



A((,fl 3 (i))<(l + ^-) 2 <C. 



K(R) <£Ml^ci^l. (7,0) 
Now we use equations (|6.6p and (|6.5p to derive 

9t(AKR) = -\rK- A + RKAP + ^(1 + ^ + A) + ^ , (7.11) 

where we used (|6.ip and (|6.2|) to substitute the second order derivatives of 
^4 and a. Since 

■ -AH, < M 



we have for sufficiently large R that 

,4 2A 2 ~ 2 

Thus we get 



Ar A% 1 
> — - H % > 



0< ^ <d t {AKR) + -RK 2 A-RKA^-^^. (7.12) 

In view of the decay estimates for A, Ar and and the bound of /3, we 
obtain 

rt a R (s,R) 

as — > as R — > oo. 



□ 



o A(s,R) 
This completes the proof of the lemma. 



8 A regularity result in maximal-isotropic coordi- 
nates 

Theorem 3 Let < e < R%. Consider a solution of the spherically 
symmetric Einstein- Vlasov system, launched by initial data in J C (R\), 
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on its maximal time interval [0,T[ of existence. If f(s,R,-,-) = 0, for 
(s,R) G [0,i[x[Q,e], then 

Ct cu+t 2 ) 

gt<Q0+9eV, (8.1) 
In particular, if f(t, R, •, •) =0 for (t, R) G [0, T[x [0, e], then T = oo. 

The last statement in the theorem is a consequence of the continuation 
criterion derived in [20]. Indeed, the requirements are that Q(t) and A(t, 0) 
are bounded. Now, since there is no matter in the domain R < e, a bound 
on A(t,0) follows from equation (|6.ip and (1T.3H . cf. the bound (|7.8p . The 
result of the theorem is not new, it is included in [20j . but as was mentioned 
in the introduction the method is different and the terms that need to be 
estimated, in order to obtain global existence without assuming a lower 
bound of R, are more regular in our approach. 

Proof: The method of proof is to a large extent analogous to the proof in 
Schwarzschild coordinates. Consider the quantities G = E(t, 1Z, W) + W > 
and H = E(t, TZ, W)- W > 0. Along a characteristic (R(t),W(t), L) we have 
by (IBTTHl) and (|6X7j) 

dG r a R aKW^ LaK aL 1 | A R ^ 

1^ IP /12-7T.2 + AZ<r>2 zrA-r. + A )i 



and 



ds 1 A E 1 2EA 2 1Z 2 A 3 1Z 2 E 1Z A 



dH r a R aKW, LaK aL 1 A R 



ds L A E 1 2EA 2 TZ 2 A 3 1Z 2 E y 1Z A 
First we consider the quantity H. Using equation (|7.1|) we conclude that 

, LaK aL 1 A R ,, C^jL C 
1 -)\ < — - — < — 

' Vt> ~ A > — T>1 — ,2 



1 2EA 3 1Z 2 A 3 1Z 2 E 1Z A ] ~ TZ 2 ~ e 2 ' 
Thus 

H{t) < H(0)efo ^-*W*s + C n ^^L dsdT _ {8A) 

Jo e 

Let us consider the first of the integrals in the exponent, the second is 
analogous. Since 

dlZ _ a W 
~ds~~ { A^~ Ph 
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we can write the integral as a curve integral 

J —(KA)(t, R)dR + - (PKA)(t, R))dt, 

where 7 is the curve (s,lZ(s)), < s < t. Let T denote the closed curve 
r := 7 + C t + Cqo + Co, oriented clockwise, where 

C t = {(t,r):R(t)<r<R OQ }, (8.5) 
Coo = {(5,^00) :t>s>0}, (8.6) 
C = {(0,r):i?oo>r>i?(0)}. (8.7) 

Here R^ > R2 + 4t, so that / = when r > i?oo- By applying Green's 
formula in the plane we get 

£-(KA)(t,R)dR+(^^ - (J3KA)(t,R))dt 

= J J dt (—(KA)(t, R)) - d R ~ (PKA)(t, R)j dtdR =: In- 

Now we use equation (16. 6p to substitute for dtK above. We obtain 

f f ,2A RR 2{A R ) 2 2A R , . , 

In= aA{-^- + ^+A^{p + p-2p T ))+aK 2 AdRdt. 

(8.8) 

From (16. ip we derive 

A RR (A R ) 2 A R 3 2 

We thus get 

U = J j n aA (—-R§-^- A7:{p - p + 2pT) y (8 - 9) 
Remark 2: The quasi local mass m is given by 

m = -(1 - |Vrl 2 ). 

Here r is the area radius, and Vr is the gradient of r. In our case r = AR, 
which implies that 

R 3 A 3 K 2 2a R 3 (A r ) 2 , 
m= R 2 A R ^f-. (8.10) 
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It thus follows that In can be written as 
r r 2m 

lQ= J J n aA (j3 R l- 47T (p-P +2 P T ')) dRdL ( 8 - U ) 

This is thus identical to the structure of the corresponding term in 
Schwarzschild coordinates, cf. Remark 1. Here we will however stick to 
the form (|8.9p for the estimates. 

To summarize, we have in view of (|8,4p obtained the estimate 

H(t) < H(0)e In - Ic t- Ic °°- 10 , (8.12) 



where 



rRoo 

I Ct = - (KA)(t,R)dR, (8.13) 
Jn(t) 



rRoo 

I Co = / (KA){0,R)dR, (8.14) 
Jn(o) 



and 



ICoo = f J( S '*i - (f3KA)(s, Roo ))ds. (8.15) 
Jo A(s,R 00 ) 

We now invoke the a priori bounds derived in the previous section. Since 
7^( s ) > e, < s < t, we have in view of Lemma EJ 

f f 2m 2t 
Iq< aA—r—rdRdt < — . 



A 3 R 3 



Using the bounds ([71]) . ([75]) and (I7TD]) we obtain 

C(l + t 2 ) 

I c t < • 

e 

The integral Ic x vanishes in view of Lemma [Hand Ic depends only on the 
initial data. Since the second term in (|8.4p can be estimated in the same 
way we obtain 

ff(t )<p(o ) + ^ a ? a . 

The estimate for G is analogous. Thus 

Ct ca+t 2 ) 
Gt<G0+9eV, 
e z 

and we get 

Ct ca+t 2 ) 
e z 

This completes the proof of Theorem [3] 

□ 
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